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We derive bulk asymptotics of skew-orthogonal polynomials (sop) ■Km \ /3 = 1, 4, defined w.r.t. 
the weight exp{-2NV (x)) , Vix) = gx'^/4 + tx^ /2, g > and t < 0. We assume that as m, N ^ oo 
there exists an e > 0, such that e < (m/N) < X^r — e, where Acr is the critical value which separates 
sop with two cuts from those with one cut. Simultaneously we derive asymptotics for the recursive 
coefficients of skew-orthogonal polynomials. The proof is based on obtaining a finite term recursion 
relation between sop and orthogonal polynomials (op) and using asymptotic results of op derived 
in Finally, we apply these asymptotic results of sop and their recursion coefficients in the 
generalized Christoffel-Darboux formula (GCD) to obtain level densities and sine-kernels in the 
bulk of the spectrum for orthogonal and symplectic ensembles of random matrices. 

PACS numbers: 02.30.Gp, 05.45.Mt 



1. INTRODUCTION 

Skew-orthogonalpolynomials are useful in the study of orthogonal (/? = 1) and symplectic (/3 — 4) ensembles of 
random matrices @, |3, 0, H, H, 0, H, H, [13, [HI [13, [ll|- this paper, we derive asymptotics of skew orthogonal 
functions (j)m\x) and ipm^x) and their recursion coefficients [1, [^, defined w.r.t. the weight 

w{x)^eM-^NVix)), V{x) = ^ + *-^, g>0, t < 0. (1.1) 

Here, 2N is a large parameter which, in the context of random matrix theory, is the size of the matrices. 
We define skew-orthogonal functions: 



1 " 
cl^i^Hx) = ^=7Ti^\x)eM-NV{x)), 7^i«(x)=5]c^^)x^ /3 = 1,4, (1.2) 

(1.3) 



i^^Hx) := </>f )(:,), ^(i)(a;) := / 0W(y)e(x - y)dy, e(r) = n € N, 
where f/i'''' are normalization constants. They satisfy skew-orthonormality relations: 

(p^^\x)^l,^^\x)dx = Z„,„, Jji'\x)^'£>{x)dx = Z^(^\l^+..., n,m G N. (1.4) 

Using these polynomials, we study the corresponding random matrix model: 

Pp.N{H)dH ^ -^exp[-[2W(ff)]]di/, /3 = 1,4, (1.5) 
where the matrix function V{H) is a double well quartic polynomial of H and 

» 2N-1 

ZpN-.^ / exp[-[2Tru(i7)]]di? = (2iV)! [] gf\ (1.6) 
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Here, M^fj^ is a set of all 2N x 2N real symmetric (/? = 1) and quaternion real self dual (/3 = 4) matrices. dH is the 
standard Haar measure. 

To study statistical properties of such matrix models, we need to study certain kernel functions 

2N-1 

For quartic potential, /3 = 1, this is given by § 

(X - y)S^2^^{x,y) = R2N-4,2Nbp2N~3i^)^2Niy) - [x y)] + R2N^2,2N+2bp2N~lix)i^2N+2iy) ^ {x y)] 

- R2N -3,2N+li''^2N -4ix)^2N+liy) — [x ^ y)] — l,27V+3 ['/'27V-2 (^) V'27V+3 (2^) — {x 2/)] 
+ R2N-2,2Nbp2N-lix)'4'2Niy) ^ (x ^ y)] - yP2N ~3,2n[^2N -4iy)'^2N (x) ~ {x y)] 

- [R2N^l,2N+lbl^2N-2ix)i^2N+liy) - {x ^ v)] - yP2N ~l,2N+2bp2N -2ix)'ip2N +2iy) - {x ^ 2/)]] 

- yP2N~2,2N+l\-'l'2N~li^)'J^2N+liy) - (x ^ V)] + yPjW- 1.27V [V'2jv-2 (2;)V'2W (j/) - (x v)] (1-7) 

where the recursion coefficients pj^^ and r'^^I are defined as: 



(i):= j;P,!;;t);'>(i); ^'>(i);=j;Rg#m(i), (L8) 



For /? = 4, ipj^\x) and (^^"'^■'(a;) are replaced by ip^^^x) and -0j'*-'(x) respectively Q. 

To study asymptotic behavior of these skew-orthogonal functions and their recursion coefficients Pj'^^ and pj^\ we 

expand (j)m\x) and il)m\x), m > 1 in a suitable basis such that their derivatives exist and 4)m\x) are polynomials 
of order m. We obtain recursion relations between sop and the corresponding op. Solving the recursion relations and 
using asymptotic properties of op [1], we derive asymptotics of sop. Finally, we apply them in the GCD formula (jl.7p 
to study the corresponding matrix models. 

A. SOP and Orthogonal Ensemble 

For sop with quartic weight, we expand Tprn (x) = 7n™Vm-3(2;) + ^k"^4''k\x) such that [V'nP]'(a;) = <^rn (x) 

is a polynomial of order m. Using skew-orthonormality (jl.4p . we get recursion relations: 

V'2m+l(2;) = 72miVV2m-2(2;) +72m-VV2r«-l(2;), m>l, (1.9) 

^2]lix) = lilTU'il-six) + 7iZ%^2l-2ix) + l2Z^4^^^ m>l, (1.10) 

AD . ^ (2m+l) rp(2) ,(2) , . p(2) A2) . J (2r„+l),(l) , . , ^ 

1'2m+l\X) — l2m+l ^27n-2,2m+l'P2m+l\X ) ^ ■ ■ ■ -^2m-2,2m-5V2m-5\X ) ^ l2m-l V2m-1\X)^ \^-^^) 

AD , \ _ {■i^)\p(.'2) ^(2)^2.^, p(2) J,(2) fa-^l + ^(2m) , (1) . . (2m) (1) , . . ^ 

Y2in\X) — 72m ^2m-S,,2m.V2m\X ) ^ ■ ■ ■ ^2m-3,2m-&V2m~&\X ) ^ l2m-2Y2m-2\X ) ^ l2m-4'P2m-4\X ) i {^■^^) 



Here 



cj)f\x) = i^i±=liexp(-W(a;)), f ^f\x)4^\x)dx = 5,-^. (1.13) 



are orthogonal polynomials. For even weight, with Rm = (hm/hm-i), Ro = 0, op satisfy [l|: 

, i+3 

x<pf\x)^,/Rj^^^^;^,{x) + ^<pf},{x); ■^^f\x)= ^ P,'lVi'^(x) (1.14) 



= -Ngy^Rj+iRj+2Rj+3, Pjj+i = "■"TTT' ^ij-i = ^l?-3 = Ng./R~^R~^j. 
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Furthermore, using (/'2m-i(0) = in (|1.14p . we get: 

4'^(o) = (-If 



R 



2m-l • ■ 



R 



R 



2m ■ 



.R 



Using skew-orthogonality relations ((</'2m-2(2^)' V'2m+i(^)) 



= 0), ((4i^_3(x),42(^)) = 



= 0) in (1111) - we get: 

(2m+l) (2m-2) p(2) 



(2m+l) 
72m- 1 

(2m) 
72m-2 



-to. 



(2m) 
72m-4 



(2m) (2m-3) p(2) 
y2m /2m-3 2m-6.2m-3 



(2m) 
72m 



(2m-l) p(2) (2m-l) (2m-3) p(2) 

72m-l ^2m-4,2m-3 + 72m-3 72m-3 ^2m-6,2m-3 



respectively. We choose 72^V^ = 1 and 72m-V^ = - V^2m-2/i?2m-3- 7^™ 2^'' 72^^2 and 



(1.15) 
0) and 

(1.16) 



from (|1.14p and (|1.16p . Here, we note that choice of 72^^^^'* and ^2ni-^i' is such that we can use properties of 
orthogonal polynomials to obtain V'2m+i(^) froni |L9]). However, this makes the sop non-monic. 
Using (fTM]) and ([06]) . we solve dM]). V'2m+i(a;) while 

'?^2m('^) is obtained using skew-orthonormalization (|1.4|): 



2m -2 
(2m+l) 



I'im-i '^an be calculated 



4m+i(2;) = 



m — 1 



^m-il^;) 



^a) = ^^LiM, ^ > 1. 



(1.17) 



2m -1 



'/'2m+i(^) and "02711 (^) <^an be obtained by taking derivative of '02m+i(^) and using skew-orthonormality respectively. 
Now, we derive asymptotics for sop using results for op [l[. In the limit m, N 00, we assume e > 0, such that 

e < A(= m/N) < Xcri— t^/'^g)- This ensures that (j)'£\x) and Vm^(cc) are concentrated on two intervals [—3:2,— a;i] 
and [xi,X2] and exponentially small outside. Using results from [ij for large m, 



^'^\x) = ?^^[cos(/„,(0)) + 0{N-% U9) = 



m+ l/2\ /sin(26') 



-0)-(-l)™| + J, (1.18) 



R2m+1, R2m — R, L + 0{N 



-t ± ^t^ - 4Ag 
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with A = gRL, t + g{R + L) = 0. 



(1.19) 



Thus we have for xi,2 = \J {—t =F 2\/Xg)/g, in the range a;i-|-i5<a;<a;2 — (5, (5>0 and m > 1, 

C'2m- 



42+l(^)^2C2m- 

where 



R 



COs(/2m-l(e))+0( ^ 



NX 



Lx 



sm 



sin(/2m-i(0))+O( ^ 



xsinO 



sin(/2m-i(f?)) + 



4'iix)^2C2rn-l 



RsinO 



C0s(/2m-l(e))+0 



q = COS t 



gx^ + 1 



cosx 



2^/Yg-tq 



2 VA^ 2 VA^g - t 

Here, we note, that for a given a;, 9 varies with to. We have chosen 



A--, C, 



1/4, 



20FVA/ 
2m -1 such that 



{\ + 0{N-^)). (1.20) 



(1) 

2m+l±k 



(x) ~ 2C2m-l 



i? 



a; sm 

and so on for i^iH^^ix), 4m+i±fe(a;), 4m±fc(^) 

For the recursion coefficients (|1.8p . we simply read off from (|1.9p - (|1.12p . We get for large m: 



[cos (/2m-i(0) T (fce)/2) + 0(l/iV)] , m » fc. 



(1.21) 



p(i) 

-^2m,2m-|-3 




-^2m+l,2m+4 ^ ; 



(1) 



■2m,2m+4 



R. 



(1) 

•2m+l,2m+5 



-A^A; 



2m,2m+l 
(1) 



5^' 



p(i) _ 

-^2m+l,2m+2 ^ 



rC. 



2m,2m+2 



(1) 

'2m+l,2m+3 



NtVRL, 



(1.22) 
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where a 



To calculate S'^]^(a;,y), we use (IL^ and in ([LTD. The first four terms (modulo 0{N~'^)) in ([LTD give 

— (2cos26'sin(a2Ar-i))/(7rsin^ 6*), the next four terms give (2 cos^ 6* sin(Q;2Ar-i))/(7r sin^ 0) while the last two terms 
give (sin(a27V-i))/7r, where aj — A9{dfj{6))/d9). Combining all, we get for x = y + Ay 



{x -y)S^^^^{x,y) 



sin(Q;27v-i) 



0{N-^) 



*s'2jv('^i y) 



sm 



Ay2N^y^l - 



irAy 



+ 0{N-'^). (1.23) 



To obtain the level-density, we take the limit: 



4m 2^5«(x,y) = ^Si%,y) = ^\y\V^+0{N-') = J^y^5 - (^)' + O(iV-), (1.24) 



while S2j!i{x,y)/ S2]y{y,y) gives the sine-kernel: 



S2Ni^^y) _ sin Trr 



AySi'^{y,y). 



(1.25) 



2. SOP AND SYMPLECIC ENSEMBLE 

For sop with quartic weight, we expand (f)m\x) = 7n™Vm''(a;) + 7fc'"Vi^H2;). Using skew-orthonormality 

(II. 4p . we get recursion relations: 



6^^^ (x) - ^(2m-H),(2) , X (2m-Hl) ,(4) . x 

02^(2;) = 7£rV2m(a;) +72^^2^-2(2;) +72n"^44m-4(2;) 



^(4) / X _ (2m+l) 
4m (a;) = 72^^ 



(2) 



'/'2m+4(2;) + 



2m+l,2m+4V^2m+4 
(2) Ai^) 



(2) 



2m+l,2m-2V 



'/•2m-2(a;) +72^W2m-l(a;) 



(2.1) 
(2.2) 
(2.3) 



2m,2m+3T2m+3\ 



4'i+3(-^) + • • ■ ^i^,2™-34'i-3(^) + 7L"J24'i-2(^) + 7t"J44'i-4(^), (2-4) 



where (x) are op defined in (|1.13p and (|1.14p . 

Using skew-orthogonality relations ((-02m-2(2;), 4m+i(^)) = 0)' ((V'2rL3(2;), 4m(2;)) 



((4m-i(^),4™(^)) ^0)mm-m, we get: 



(2m+l) _ „ (2m+l) (2m-2)p(2) 
/2m-l ^ ^ /2m+l f2m-2 ^2m-2,2m+l i 



(2m) _ (2m) (2m-3) p(2) 

72m-4 ~ ^72m 72m-3 ^2m-3,2m 



(2m) 
72m-2 



„ (2m) 
^72m 



■ (2m-l)p(2) (2m-l) (2m-3)p(2) 

72m-l ^2m-l,2m ~r 72m-3 72m-3 -'2m-3,2m 



0) and 



(2.5) 



v/-R2m+i/i?2m. 72m- 2^^ : 72m-2 ^nd 72^1^4 Can bc Calculated 



respectively. We choose 73^^^' = l/\/2 and 73^^^ 

from (|1.14p and (|2.5p . Choice of 72^V'' ^-^^d 72m^V'' such that we can use properties of orthogonal polynomials to 
solve (j2.ip . However, this makes the sop non-monic. 

;,(4) 



Using (|1.14p . (|1.15p and (|2.5p we solve (|2.ip . V'2m(^) obtained using skew-orthonormalization (|1.4p : 



4tUr(-) - 



a; 



72 



411+2(2;) - 4m+2(0) CXp[-iVl-(.T)] , 4'J^ix) - 



a;4m+2(2;) 



•\/2i?2m+2 

(x) respec 

With m, N 00 and in the range xi + S < x < X2 — S, S > 0, and neglecting the m independent term, we get 



^2m(^) ^n^i V'2m+i(^) "^^^ derived by integrating and differentiating ip^mix) and 02m+i(^) respectively 



(2.6) 



4'^+i(^)^c' 



2L 



x sm 



[cos [/2m+2(e)]+0(l/A^)] 



4'^(^) 



C 



V'2m+l(2;) 



4A^CA' 



/2i? 



iVA' V 2sin3 6l 

[sin [/2m+2(e)] + 0(l/iV)] , 4m (^) " 



[sin [/2m+2(e)]+0(l/A^)], 



2a;3 



L sin Q 



[cos [/2m+2(^?)]+0(l/iV)] 



(2.7) 
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where fm{d), Cm = C, x, A and 9 are defined in (|1.18p and (|1.20p . Here, we note, that for a given x, 9 varies with m. 
We have chosen 9 = 6*2^+2 such that 



2L 
xsin^ 



cos 



0{l/N) 



» fc, 



(2.8) 



and so on for ^'^^Jx), V'sm+iiJs;), V-atlifcC^^)- 

For the recursion coefficients (jl.Sp . we simply read off from (|2.1|1 - (|2.4|1 and use (|1.14p . For large m we have: 



(4) 



2m, 2m +3 



R 



IR 

(4) 

2m,2m+4 



-^2m+l,2m+4 ^ " -'^J 



p(4) p(4) _Q 

^2m,2m+l ^ > ^2m+l,2m+2 ^ 



(1) 



■2m+l,2m+5 



5 = -NX, 



^2m,2m+2 ^ ^2m+l,2rn+3 ^ JVI,vn,lj. 



(2.9) 



To calculate S';^^(a;,y), we use ^J\, jM]) and dm in dLT]) for /3 = 4. The first four terms (modulo 0(7V"i)) in 
p.7p for /? = 4 give — (cos20sin(a2Af+2))/(7''sin^ 0), the next four terms give (cos^ sin(a2Ar+2))/(7'' sin'^ 0) while the 
last two terms give (sin(a2Ar+2 ))/"'• Combining all, we get for x = y + Ay, 



(,-.)^(^^(x,,) = -5^^%^ + 0(A^-) 
For level-density, we take the limit 



sm 



2nAy 



O(iV-i) (2.10) 



1 



Ay^O N 

such that we get the "universal" sine-kernel 



V5, 

TT 



TT 



5^^(a:,y) _ sin 27rr 



2iTr 



^ySi%,y). 



^ gy^ + t 



+ 0{N-^) 



(2.11) 



(2.12) 



3. CONCLUSION 



The key achievement of this paper is the derivation of Eqs. (|1.17p and (12. 6p . This enables us to use asymptotic 
results of op [l| to derive bulk asymptotics of sop with quartic weight [IJj, J^, JJ,] . Simultaneously (|1.9p - (|1.12p and 
(I2.ip - (|2.4p gives us the recursion coefficients Pj^j} , ^j^k- These results for sop are applied in the GCD formula [l] to 
study quartic orthogonal and quartic symplectic ensembles of random matrices in the bulk. We note that asymptotics 
for sop away from the bulk can be trivially obtained from (|1.17p and p.6p since the corresponding results for op are 
already known. 
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